Imaging of single infrared, optical, and ultraviolet photons using distributed tunnel 
junction readout on superconducting absorbers 
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Single-photon imaging spectrometers of high quantum efficiency in the infrared to ultraviolet 
wavelength range, with good timing resolution and with a vanishing dark count rate are on top 
of the wish list in earth-bound astronomy, material and medical sciences, or quantum information 
technologies. We review and present improved operation of a cryogenic detector system potentially 
offering all these qualities. It is based on a superconducting absorber strip read out with super- 
conducting tunnel junctions. The detector performance is discussed in terms of responsivity, noise 
properties, energy and position resolution. Dynamic processes involved in the signal creation and de- 
tection are investigated for a basic understanding of the physics, and for possible application-specific 
modifications of device characteristics. 
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I. INTRODUCTION 

Superconducting tunnel junction (STJ) detectors are 
among the most advanced cryogenic sensorsi with intrin- 
sic spectroscopic resolution and high detection efficiency 
over a broad energy range. Essential advantages of low- 
temperature detectors in general are the small energy 
gap in superconductors compared to standard semicon- 
ductors (typically three orders of magnitude lower), de- 
coupling of the electron and phonon systems in normal 
metals below about 50 mK, and a vanishing lattice heat 
capacity. 

The lifetime of noncquilibrium quasiparticles gener- 
ated in superconductors due to local energy deposition 
in excess of a purely thermal distribution becomes very 
large at operating temperatures well below the transition 
temperature. Those quasiparticles can be efficiently de- 
tected as a tunnel current in an STJ if tunneling rates 
are high compared to recombination and loss processes. 
In the tunneling process, a quasiparticle can tunnel back 
and contribute many times to the total charge transfer^ 
This internal amplification process is a unique feature of 
STJs which significantly improves the sensitivity of the 
detectors, allowing us to observe single sub-eV energy 
quanta. 

As a next generation detector, STJs are very promis- 
ing candidates e.g. for astronomical observations In 
addition to their direct spectroscopic response and high 
speed (^ 10^ counts/s), they can be operated with es- 
sentially no dark count rate, which is another enormous 
advantage compared to CCD imaging systems, making 
them the favored choice for detection of faint objects. 
However, to satisfy the desire for detectors with imaging 
capabilities as offered by CCD cameras, a competitive 
STJ based multipixel camera inevitably faces the com- 
plexity of large channel-number readout, i.e. the problem 
of transferring the small charge-signals from the cryo- 
genic environment to external room-temperature elec- 
tronics. One possible solution to reduce the number of 
readout lines is by using distributed readout imaging de- 



vices (DROIDs), where the quasiparticle distribution in 
a large superconducting absorber is detected with several 
STJs at different positions. This type of detector was in- 
vestigated with x-rays^iSiLSiSiifliii as well as with optical 
photonspi2iiiili^ 

In this paper we present measurements with quasi- 
onedimensional DROIDs irradiated with energy quanta 
covering the ultraviolet (UV) and the entire optical 
range, and including an extension to near infrared (IR) 
energies below 1 eV. After a description of device prop- 
erties, experimental conditions and measurement proce- 
dure in Sec. the linearity of detector response, the 
noise performance and the resolution in energy and po- 
sition are discussed in Sec. IIII Al Time-dependent pro- 
cesses in the detector involving quasiparticle diffusion, 
loss and trapping rates, are considered in more detail in 
Sec. IIII Bl The paper concludes in Sec. IIVI 



II. EXPERIMENT 

The devices were fabricated by sputter deposition on 
a sapphire substrate. The Ta absorber was deposited 
first, with an area of Lxa x = 135 x 31.5 /im^ and a 
thickness of 100 nm. The square shaped Ta-Al junctions 
with a side length Lj = 25 /xm and an Al layer thick- 
ness of 38 nm were fabricated at each end on top of the 
absorber. Hence, the length of the bare absorber was 
^b,Ta = 80 /im. The residual resistivity ratio of the ab- 
sorber was RRR = Pht/Plt ~ 25, where p^r^ and p^T 
the (normal) resistivities at room temperature and low 
temperature, respectively. 

The cryogenic detectors were operated at 310 mK and 
the supercurrent was suppressed by application of a mag- 
netic field parallel to the tunnel barrier. With two 
independent charge-sensitive preamplifiers operated at 
room temperature the devices were voltage biased at 
Vf, — 80 /iV where a thermal current below 50 pA was 
measured. At this bias voltage the current / was found 
to depend on temperature T according to the empirical 
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expression / = 0.345 exp(T/26.4 mK) fA. Although the 
minimum thermal current was observed at Vf, « 200 /xV, 
the responsivity of the devices showed a maximum at 
lower Vb, in the range of negative differential resistance of 
the current-voltage characteristics. We introduced ohmic 
resistors of 5 kil in the signal lines (in series) close to the 
devices, which acted together with the junction capaci- 
tances as efficient RC filters against high-frequency noise 
on the 1.2 m cables, as well as damping resistors against 
resonances of the junction-cable system corresponding to 
an LRC circuit with intrinsically low damping constant. 
Compared to the superconducting device resistance on 
the order of 1 Mil the influence of the filters on the de- 
tector response was negligible. 

A variety of pulsed light-emitting diodes (LEDs) was 
used as light sources, with wavelengths ranging from IR 
at 1550 nm (0.801 eV) to UV at 370 nm (3.355 eV). The 
specified spectral widths of the LEDs were typically a 
few percent. The photons were coupled to the detectors 
via optical fibres and through the sapphire substrate. 
Back-illumination through the substrate was found to 
be essential for optimum detector performance, probably 
thanks to a clean, oxide-free interface, whereas the native 
TaOx on top of the absorber appears to be less transpar- 
ent. The light intensity (pulse duration) was adjusted to 
observe on the order of one photon per pulse. Packets 
of (incoherent) multiple photons were synchronized to a 
maximum time spread of less than 200 ns. 

Detection efficiency of Ta absorbers is about 60% for 
optical and UV photons, but it drops drastically below 
red photon energies ^ 2 eViJ^^ Radiation shorter than 
about 200 nm is expected to be cut off by the sapphire 
substrate. Unfortunately, we had no reference system 
for calibrating the absolute photon absorption probabil- 
ity. 

The electronic signals from the preamplifiers had a rise- 
time of about 1 fxs (comparable to the tunneling time and 
to the ~ 1 MHz bandwidth of the readout electronics) 
and a decay on the order of 30 fj,s. The preamplifier out- 
puts were band-pass filtered with a 10 /iS time constant 
and digitized for offline analysis. 

Due to noticeable dependence of the device respon- 
sivity on Vf, and on electronics settings the system was 
calibrated against an electronic pulser injecting a defined 
number of charges. The charge noise of the readout elec- 
tronics with the junctions connected (but not irradiated) 
was determined from pulser distribution spread or from 
rms noise measurements to be q,, « 2760 e for both chan- 
nels combined, i.e. about 1950 e for each individual chan- 
nel (at Vt = SOfiY). 

Figure n shows results from measurements with UV as 
well as with IR photons. The charges Qi and Q2 refer 
to the separate outputs from the two channels. The left- 
hand plots display the total charge = Q1 + Q2 versus 
the normalized (onedimensional) position of the photon 
interaction xq — (Qi — Q2)/Qt.- Slight asymmetries of 
the data with respect to xq are due to the difficulty of 
perfect adjustment of operating point for the two chan- 



^ 225 
O? 200 
175 
150 
125 
100 
75 
50 
25 


80 

O) 70 
60 
50 
40 
30 
20 
10 








: (a) 

































(b) 




























1 






















noise 







-0.5 




-0.5 



0.5 100 

Xg events/bin (1/ke) 



FIG. 1; (a) Total charge output Qe versus normalized posi- 
tion for the DROID irradiated with 3.355 eV (370 nm) UV 
photons, (b) Projection of the charge spectrum (histogram) 
of the data in (a) for |a;o| < 0.05. The labels to the right 
of the peaks refer to the number of simultaneous photons, 
noise-triggered events are labelled with 'noise'. (c),(d) Data 
as presented in (a),(b), but with the detector irradiated with 
0.801 eV (1550 nm) IR photons. 



nels. The shape of the event distribution is determined 
(in first order) by quasiparticle diffusion, loss processes, 
local energy gap, and the efflciencics of trapping and tun- 
nefing (see Sec lIIIBll . 

Spectral response and resolution of the experimen- 
tal data are determined for events with |a;o| < 0.05, a 
projection of which is shown in the right-hand plots of 
Fig. ^ The number of simultaneously absorbed energy 
quanta is given to the right of the corresponding peaks, 
whereas 'noise' denotes events triggered on photonless 
fluctuations. A clear spectral separation of single-photon 
events from multiple photons and from noise triggers is 
observed, down to the lowest measured IR energy quanta 
of 0.801 eV 



III. DISCUSSION 

As a first remark, we wish to comment on the events 
observed in Fig. ^a) in the range (0.3 ^ ^ 0.6) and 
(55 ke < ^ 75 ke) . They originate from single-photon 
energy depositions in the ground lead connecting to the 
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absorber. This lead is attached asymmetrically to the 
detector, very close to junction 1. Therefore, a fraction 
of the quasiparticles generated in the ground electrode 
leak into that STJ and contribute a measurable signal. 
This effect exemplifies that one can in principle design 
any absorber shape appropriate for a specific imaging 
purpose. 

The expected number A^o of excited quasiparticles gen- 
erated upon deposition of energy E is Nq = £'/!.? A 
where A is the energy gap of the superconducting ab- 
sorber and 1.7 A is the effective energy to create one 
quasiparticleii2iSSi2i Note that the energy conversion fac- 
tor 1.7 is appropriate for Ta, but in general it is material 
dependent iSl The intrinsic energy fluctuations Si of STJs 
can be described in first order by^^i 



= v/l.7A(F+l + (n)-i), (1) 

where F « 0.2 is the Fano facto]»l^^, and (n) is the aver- 
age number each quasiparticle contributes to the charge 
signal due to backtunnelingi^ The energy- resolving power 
R is conventionally described by i? = E/2.355ei where 
the factor 2.355 scales between standard deviation and 
full width at half maximum of a normal distribution. 

A. Linear response and spectral resolution 

At temperatures well below the superconducting tran- 
sition temperature, where quasiparticle recombination 
processes are very slow, and for sufficiently low energy 
densities of the nonequilibrium quasiparticle distribution 
(i.e. for optical photons, where self-recombination is neg- 
ligible) , a linear response of the detector is predicted for 
single STJ detectors^ as well as for DROIDsi^i In or- 
der to test the linearity of our detectors we performed 
measurements with photon energies ranging from IR to 
UV. Similar to Figs, ^b) and (d) we consider only pho- 
ton events with interactions in a narrow window of the 
absorber's central region, i.e. satisfying |a;o| < 0.05. In 
addition to single-photon events we also extract the sig- 
nals due to simultaneous two-photon events where data 
of minimum interaction distance is selected. This con- 
dition corresponds to the low edge of the two-photon 
events contour in Fig. Q^a) and to the prominent peak 
at ^ 160 ke in Fig. ^b), whereas events of two photons 
interacting at largest distance (i.e. close to the junctions) 
are expected to appear at the upper contour edge. 

The results of detected total charge Qs as a function 
of photon energy are shown in Fig. [2Ja) . Solid dots and 
open circles correspond to single and two photons, re- 
spectively. The horizontal error bars refer to the spec- 
tral widths e\ of the LEDs, the vertical errors reflect the 
measured distribution spreads qm of signal amplitudes. 
A linear least-square fit to the data points yields 

r^^Q^/E^ 23843 e/eV. 

This is by a factor 28.4 times more charge output than the 
theoretically generated Nq = 840 quasiparticles per eV in 
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FIG. 2: (a) Detected total charge as a function of photon 
energy. Full dots and open circles correspond to single and 
two-photon events, respectively. The horizontal error bars are 
the specified spectral widths of the LEDs, the vertical errors 
are taken from the experimental detector signal distributions. 
The straight line is a least-square fit to the data points with 
a slope as indicated, (b) Intrinsic experimental detector noise 
(electronic noise and photon energy spread subtracted) versus 
energy. The drawn line is a theoretically predicted thermo- 
dynamic limit. The right ordinate relates the data to units 
of eV. (c) Ratio of measured to theoretical intrinsic detec- 
tor noise, (d) Estimated onedimensional position resolution, 
(e) Energy-resolving power R = Qs/2.355qi. 



a Ta absorber, where Axa = 0.7 meV. This amplification 
is attributed to the device-internal gain due to backtun- 
neling. In the case of lossless diffusion the responsivity 
would even amount to 27.9ke/eV (see Sec. IIII B|l . corre- 
sponding to an average backtunneling factor (n) ~ 33.2. 

We found no significant deviations from linearity 
within experimental errors over the entire measured en- 
ergy range including the data from two-photon events. 
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This is most remarkable because it proves not only the 
proportional regime considered in Refs. 12312^ but also 
an energy-independent backtunneling process at these 
quasiparticle densities. 

By subtracting the readout electronics noise Qe and the 
spectral widths e\ of the LEDs from the measured pulse- 
height distribution spreads qm we deduce the device- 



intrinsic charge fluctuations qi as qf 



Figure E^b) shows qi as a function of energy together 
with the theoretical fluctuations (?th = Si'H given by 
Eq. iQ}. The ratio g^/gth ~ 1-4 plotted in Fig. |3c) sug- 
gests that the performance of our detectors is close to the 
thermodynamic limit imposed by the simple model l|T]l. 
The constant ratio implies that further additive, energy- 
independent terms in Eq. may better account for the 
device-intrinsic fluctuations. ^^'^^ 

If we neglect quasiparticle losses during diffusion, 
the charge noise can be translated into an imaging 
position resolution 5x « (-^Ta/a^max)(9i/Qs) as dis- 
played in Fig. Old), where Xmax = 0.41 in our de- 
vices accounts for the reduced trapping efficiency (see 
Sec. IIII B|) . The energy-dependent number of effective 
pixels ~ LTa/2.355fe is in our case in the range of 3 ... 6 
(neglecting the real geometry where about one third of 
the absorber is modified by trapping regions). A possi- 
ble improvement in position resolution, however, can be 
obtained for longer absorbers (for the price of reduced en- 
ergy resolution due to diffusion losses) by measuring the 
time delay of the charge signals in the two channels. As- 
suming a moderate timing accuracy of 1 /xs we estimate 
from our timing measurements (roughly 6 /iS differential 
delay over Lb,Ta, see Sec. llllBTl a spatial resolution of ~ 6 
pixels, which is already comparable to the charge-noise 
limited values. 

Finally, the energy-resolving power R of our DROIDs 
is plotted in Fig. [2Ie), showing a comfortable signal-to- 
noise ratio for single-photon imaging down to near IR 
wavelengths. 



B. Quasiparticle diffusion 

After localized energy deposition in a superconductor, 
cascades of several fast processes eventually end with a 
population of excess quasiparticles. Those initial con- 
version processes are negligibly fast compared to quasi- 
particle recombination and loss times Therefore, we 
expect our measurable device dynamics to be dominated 
by diffusion and subsequent trapping and tunneling pro- 
cesses. Quasiparticle diffusion in the absorber can be 
described by the onedimensional equation 



d_ 

dt 



1 

n(x, t) - D—^n{x, t) H 

dx^ Tloss 



nix,t)=0, (2) 



where n{x, t) is the quasiparticle density, D is the diffu- 
sion constant, and t^^^ is the loss rate in the Ta absorber. 
We adopt the phenomenological model in Ref. 6 which 
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FIG. 3: Contour plots of a measurement with A — 592 nm 
photons (only single- photon events are shown). The thick 
line represents the result of a model fitted to the data points, 
where we assume onedimensional quasiparticle diffusion with 
losses and a reduced gap in the trap region of the Ta absorber 
due to proximity of Al. Photon absorption in the trap region 
is identified with events at xq — ±a;max. 



derives the final integral of collected charges Qk=i.2 at 
the two STJs 



]k{xo) = Qo- 



sinh^fea + /Scosh^^Q; 
1 + ) sinh a + 2(3 cosh a 



(3) 



where — \xo ^ (^l)''l/2 is a normalized photon po- 
sition relative to junction fc, the dimensionless parame- 
ter a — L / \/ Dt\oss measures the quasiparticle diffusion 
length A = L/a relative to the absorber length, and 
/? = (Ttiap/Tioss)^^^ compares the trapping rate t^.^^ to 
the loss rate. 

This model allows us to fit the experimental data as 
shown in Fig.|2| In our specific devices, we have to distin- 
guish between the bare Ta absorber and the proximized 
region due to trapping layers, where the effective gap 
energy relevant for quasiparticle generation is reduced. 
Those regions are separated at \xo\ w 0.3. The position 
identified with energy depositions directly in the trap re- 
gion is Xmax = 0.41 in our DROIDs. Because a accounts 
for the curvature (losses) and f3 for the given Xmax (trap- 
ping), we can fit Eq. iPJ to the data in the limited range 
|a;o| < 0.3, yielding the fit parameters a = 0.63, /3 = 0.57 
and Qo,b = 58.6 ke for the bare absorber. Applying the 
same function to the data for |a;o| > 0.3 while keeping 
the formerly found a and (3 fixed yields Qo.t = 69.1 ke 
in the trap regions. This implies that the mean gap en- 
ergy in the trap region where quasiparticle generation 
takes place is reduced by (1 — Qo,b/Qo,t) = 15.2%. The 
empirical fit to a Monte Carlo simulation^-- has found a 
reduction of about 12.5% for the same devices. The dis- 
crepancy between the results from the two approaches 
is acceptable considering the uncertainties in empirical 
parameter adjustment in the latter method. One should 
note that the energy gap of 0.848ATa = 594 /xeV in the 
trap region is not equivalent to the one at the tunnel bar- 
rier where about 450 /leV was found from current-voltage 
characteristics measurements. This is well understood in 
terms of the superconducting proximity effect 
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For the parameter 13 the simulations-'^- found j3 « 
(-Pioss.Abs/^trap)^/^ = 0.483. This value is sUghtly lower 
than our model fit parameter because the simulations^^ 
consider a trapping and a tunneling probability sep- 
arately (i.e. a nonvanishing detrapping probability), 
whereas our model averages over trapping and detrap- 
ping processes and therefore results in a lower effective 
trapping rate. 

The quasiparticle diffusion length in our absorber, i.e. 
the average length a quasiparticle travels before being 
lost by recombination or other loss channels, is 

A = — = 2f4^m, 
a 

which suggests that there is still room for a longitudinal 
extension of the absorber. In the expression above we use 
LTa (and not Lb,Ta) as the effective absorber length based 
on the following argument: A quasiparticle entering the 
trap region travels, if not tunneling before, a distance 2Lj 
until it has the first chance to detrap. The average trap- 
ping length A/3^ is on the order of 2Lj (calculated with ei- 
ther Lxa or Lb.Ta)- The mean distance in the trap region 
travelled by quasiparticles which are being trapped is 
therefore roughly Lj. Consequently, the absorber length 
relevant for our analysis is Lb,Ta + 2Lj = Lxa- Alter- 
natively, we can consider quasiparticles starting in the 
trap region with arbitrary direction and position. Those 
which are detected in the opposite STJ travel in aver- 
age a distance Lj before leaving the trap. Therefore, the 
events at |xniax| in Fig. |3| are identified with the phys- 
ical positions at ±LTa/2, and not with the edge of the 
transition from bare absorber to trap region at ±Lb,Ta/2. 

Measurements of the trigger time differences 5t,n be- 
tween the two pulse-shaped detector signals for thresh- 
olds set arbitrarily at 33% and 67% of the peak ampli- 
tudes are shown in Fig. ^ An S-shaped timing versus 
position curve, which we observed for fast signals, was 
stretched to a straight line after pulse-shaping with a 
time constant one order of magnitude larger than the 
original signal rise-time. The results of a linear fit to the 
data points in the range |xo| < 0.3 are included in the 
figure. The events corresponding to photon hits in the 
trap regions are excluded from the fit because they show 
a slight extra time delay. 

This observation is interpreted as a consequence of de- 
layed detrapping, namely the quasiparticles detected in 
the opposite STJ have a non-zero probability to propa- 
gate within the trap for some time before leaving for the 
other side. In addition, they may initially go through a 
few backtunneling cycles, which further delays the oppo- 
site charge signal by the extra tunneling times. 

We have numerically simulated the quasiparticle prop- 
agation according to Eq. ||2Jl to determine the detector- 
signal dynamics including the response of readout elec- 
tronics and pulse shaping. The normalized timing differ- 
ences 5tn between the two channels calculated for trigger 
thresholds at 33% and 67% were 0.336 and 0.351, re- 
spectively, assuming D = L = 1. Scaling the measured 
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FIG. 4: Measured time differences between the two detector 
channels versus position, for trigger thresholds at 33% (top) 
and at 67% (bottom) of the signal peak amplitudes. The 
white lines are linear fits to the data in the range |2;o| < 0.3. 
The resulting fit parameter is included in the plots. 



data with these results allows us to estimate the diffusion 
constant 

i^«-i^i|,«7.8^. 

For comparison, a theoretically evaluated diffusion con- 
stant is given bjiSS 



D 



th 



27rA 



VpTh 



is the Fermi velocity of 



where vp = 1.6 x 10^ 
electrons in the superconductor, ri,„p = 171^/ Pjjj,ne^ is 
the impurity-scattering time, rUe is the electron rest 
mass and n — 5.52 x 10^^ cm~'^ the density of conduc- 
tion electrons in Tai^ With p 



= 13 X 10-s Om and RRR 



= 0.5 X 10-**flm (or 
25, respectively) and 
at T = 310 mK we obtain Dth = 25.2cm^s~^, which 
is more than three times the experimentally determined 
value. Similar discrepancies between experiment and the- 
oretical predictions were systematically observed in other 
experiments^&iSiiiiSSi^i and remain unresolved. 

The quasiparticle loss time deduced from the measure- 
ments is found to be tioss — A^Z)^^ w 58.9 /iS, which 
is significantly longer than the tunneling time of about 
1 /is. However, the trapping process with a time constant 



''"trap — /3 Tlo 



18.9 ^s should preferably be improved 



in our devices towards faster trapping rates, in order to 
enhance Xmax — *■ 1 for better position resolution. 

Alternatively to our phenomenological discussion on 
imaging resolution in the previous section, the resolving 
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capabilities have been derived analytically^A with a pre- 
diction for 

5x « V Dtioss = A-^ . 

Inserting our experimental results (always for single- 
photon events only) yields position resolutions in the 
range Sx w 26.7 . . . 14.8 /xm, corresponding to 5 ... 9 pix- 
els for our detector. This is in fair agreement with our 
former approximations. However, we wish to empha- 
size that all these rough estimates ignore the geometrical 
properties of the real absorber, which needs to be taken 
into account for a quantitatively precise analysis of posi- 
tion resolution»H 



IV. CONCLUSIONS 

Detection of single (and simultaneous multiple) pho- 
tons with good spectral and spacial resolution was per- 
formed with ST J based DROIDs. Single-photon resolu- 
tion down to near IR energies was proven, with a per- 
fectly linear energy response over the entire UV to IR 
wavelength range. Sensitivity to the investigated photon 
energies was possible due to the backtunneling effect, de- 
livering in our case a device-intrinsic gain of about 33, 
which is a feature unique to STJs. The detectors were 
found to operate close to the thermodynamic limit im- 
posed by particle fluctuation statistics. 

The measured dynamic response of the DROIDs was 
compared to numerical modelling based on quasiparticle 



diffusion including loss and trapping processes. The pa- 
rameters extracted from experimental data are physically 
meaningful and coincide reasonably with former Monte 
Carlo simulations.^^ However, the diffusion constant was 
found to significantly disagree with theoretical predic- 
tions, similar to all preceding experiments of the same 
kind. 

Position resolution of our relatively short absorber was 
estimated to about 3 ... 6 equivalent pixels for photon 
energies in the range 0.8...3.4eV. By taking the dif- 
ferential time delays of the two channels into account, 
we predict an improvement in position resolution by at 
least a factor of two for a longer absorber (Lxa — ^ A) and 
better trapping efficiencies {rtrap ^ ^ioss)- 

High-sensitivity spectrometers with single-photon 
counting capabilities in the broad optical range are 
not only of high interest for astronomical observations. 
Single photons at the telecommunication wavelength 
of 1550 nm are currently used, e.g., in intense studies 
on quantum cryptography^^ which links information 
theory to quantum entanglement physics. Hence, STJ 
based DROIDs with the properties and potential as 
presented in this paper may as well be an interesting and 
natural choice for IR single-photon counting experiments. 
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